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ABSTRACT 

Context. Tidal interactions of galaxies in galaxy clusters have been proposed as one potential explanation of the morphology-density 
relation at low masses. Earlier studies have shown that galaxy harassment is a suitable mechanism for inducing a morphological 
transformation from low-mass late-type disk galaxies to the abundant early-type galaxies. 

Aims. The efficiency of tidal transformation is expected to depend strongly on the orbit of a galaxy within the cluster halo. The orbit 
determines both the strength of the cluster’s global tidal field and the probability of encounters with other cluster members. Here we 
aim to explore these dependencies. 

Methods. We use a combination of N-body simulation and Monte Carlo method to study the efficiency of the transformation of 
late-type galaxies by tidal interactions on different orbits in a galaxy cluster. Additionally, we investigate the effect of an inclination 
between the disk of the infalling galaxy and its orbital plane. We compare our results to observational data to assess the possible 
relevance of such transformations for the existing cluster galaxy population. 

Results. We find that galaxies that entered a cluster from the outskirts are unlikely to be significantly transformed (stellar mass loss 
< 6%). Closer to the cluster centre, tidal interactions are a more efficient mechanism (stellar mass loss up to 50%) for producing 
harassed galaxies. The inclination of the disk can reduce the mass loss significantly, yet it amplifies the thickening of the galaxy disk. 
Galaxies with smaller sizes on intermediate orbits are nearly unaffected by tidal interactions. The tidal influence on an infalling galaxy 
and the likelihood that it leads to galaxy harassment make a very stochastical process that depends on the galaxy’s specific history. 
Conclusions. We conclude that harassment is a suitable mechanism that could explain the transformation of at least a fraction of 
galaxies inside galaxy clusters. However, the transformation would have to start at an early epoch in protocluster environments and 
continue until today, in order to result in a complete morphological transformation. 

Key words. Galaxies: evolution - Galaxies: interactions - Galaxies: dwarf - Galaxies: structure - Galaxies: kinematics and dynamics 
Galaxies: clusters: general 


1. Introduction 

For bright galaxies, iDresslerl ( 1980 ) showed that the popula¬ 
tion fraction of elliptical galaxies in clusters increases with in¬ 
creasing local environmental galaxy density, while at the same 
time the fraction of spiral galaxies decreases. An analogous 
morphology-de nsity relation is observed for dwarf galaxy types 
feinggeli et al.1 Il987l) . Dwarf irregulars and low-mass spirals 
dominate outside of dense environments, while early-type dwarf 
galaxies (dE, subsuming objects termed dwarf lenticulars) are 
the most common type of galaxy in clusters. However, the ori¬ 
gin of this large population of dEs is still a topic of ongoing 
debate, partly because they were found to follow a number of 
continuous relations in structure a nd colour with bright ellipti- 
cal and/or lenticular galaxies (e.g . IC.ote et aTll2007b ICfaTn et al l 
l2010t iKormendv & BendeH 120121) . but also since their identi- 
fication a nd d etailed ana lysis are not yet possible at higher 
shifts (see iBarazza et al . (1200 9) for dEs in a z=0.165 cluster). 


* e-mail: TL@x-astro.net 


Clues to the origin of dEs may be given by the prese n ce o f 
weak disk features in a significant fraction of the m (IJerien et al.l 
l2000t iDe Riicke et al.l 120031: iLisker et ai1l2006bl) . the presence 
of residual star fo rmation activity in dEs outside of the dens¬ 
est environmen ts (ILisker et al.ll2006at iTullv & TrenthamlfeOOSl: 
iPak et al.ll2014l) and the finding that these have significantly flat¬ 
ter shapes than d Es w ith low relative velocities in the Virgo 
cluster core (iLisker et al .1120091) . Taken together, this may be re¬ 
garded as supporting the hypothesis that many dEs were formed 
from disk galaxies that were affected and eventually transformed 
by a dense environment (for f urther observational aspects re- 


lated to this hypothesis, see also KormendvO 

I985t Binggeli et al. 

1990 

; De Riicke et al.| 120051: Beasley et al. 

2006;; Smith et al. 

2012 

: Rvs et al. 2014l). In this scenario, dEs with disk features 


dEdis) could be galaxies in an incomplete tra nsformation phase 
iGraham et al.ll2003l: iMastropietro et al.ll2005l) . 


Different processes have been discussed to explain the mor¬ 
phological transformation. In general, they can be divided into 
two main types based on the underlying physics: hydrodynami- 
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cal scenarios and tidal scenarios. In reality, a mix of both types 
will act on an infalling galaxy, but it is still an open and rel¬ 
evant question how strong the impact is for which kind of in¬ 
teraction. One of the m ost p opular tidal scenarios is galaxy ha¬ 
ras sment[Moore]etaD (1 9961) . It is characterised as the combined 
effect of tidal interactions between an infalling galaxy and the 
global cluster potential and of the tidal interaction with dose 
high speed encounters with other cluster members. [Moore et "all 
( 1996 ) used a set of N-body simulations to show that these tidal 
interactions are suitable for causing a morphological transfor¬ 
mation of spiral galaxies. This outcome is generally supported 
by studies of various authors, even though debating the actual 
efficiency dMoore et al.l Il998l Il999l: IMastropietro et al.l 120051: 

lAguerri & Gonzalez-Garciall2009l: ISmith et al.ll2Q10h. 

The efficiency of the harassment scenario is expected to de¬ 
pend on the orbit of an infalling galaxy within a galaxy cluster 
(ISmith et all 1201 Ol) . The deeper an orbit is placed inside a clus¬ 
ter, the stronger the tidal force of the background mass of the 
cluster and the higher the probability of c l ose en counters with 
other cluster members. iMastropietro et al . (12005!) show in their 
simulations of bulgeless disk galaxies with masses of 7 • 10 10 M o 
that combined tidal interactions are efficient at harassing galax¬ 
ies in the innermost regions of a cluster on orbits with pericen- 
tres very close to the cluster centre (between 40kpc and lOOkpc) 
and less efficient on orbits in the outer regions of a cluster. 
[Smith et al . (12010 ) show that late-type dwarf irregular galaxies 
with masses of 10 10 M o that enter a cluster from its virial radius 
are only slightly affected and not harassed, but they find that 
dwarfs on orbits very close to the cluster centre comparable to 
those of iMastropietro et ah (200 5) could be harassed efficiently. 
Both studies have shown what happens to late-type galaxies of 
different mass ranges on orbits very close to the cluster centre 
and far away from it. But especially on more moderate orbits, 
the question is still open whether galaxies could be efficiently 
harassed or not. With the study presented here, we intend to add 
a piece to that puzzle. 

Complementary to the orbital parameters, iMoore et"afl 
( 1999 ) show that the properties of the progenitor galaxies are 
relevant for the efficiency of transformation by tidal interactions. 
They show that compact high surface brightness galaxies (HSB) 
are hard to affect, and more extended low surface brightness 
galaxies (LSB) can be affected efficiently. A comparative result 
was shown by [Chang etaD (12013 ) by probing the efficiency of 
tidal stripping for different morphologies of satellite galaxies in 
host galaxy of Milky Way-like mass. In their simulations of disk 
galaxies orbiting inside dark matter haloes with group masses, 
[Villalobos et al.1 ( 2012 ) have additionally found that the incli¬ 
nation between the orbital plane and the disk of the infalling 
galaxy is an important parameter for the effect of tidal stripping. 
However, tidal stripping is only the most extreme effect of tidal 
interactions, and one has to expect that an inclination of the in¬ 
falling galaxy may affect the reshaping of the galaxy. 

To sum up, the parameter space for the efficiency of the 
tidally induced transformation in a galaxy cluster spans the or¬ 
bital parameters, the type of progenitor, and the orientation of 
the galactic disk. Furthermore, the properties of the cluster itself, 
such as its mass, its galaxy density, and the velocity dispersion 
of its galaxies, are expected to be relevant. Some parameters are 
partly degenerated; e.g., a variation in the orbit and the cluster 
properties could end up with a similar tidal strength. To deter¬ 
mine the contribution of harassed galaxies to the population of 
dEs in clusters, it is necessary to understand the efficiency of 
harassment depending on these parameters. 


In this study, we investigate the efficiency of tidally induced 
transformations in a Virgo-like cluster by varying the orbital pa¬ 
rameters. We adopt three different orbits. On the innermost or¬ 
bit, we use three different orientations of the disk with respect to 
the orbital plane. Furthermore, on the innermost orbit, we probe 
how the efficiency changes if the progenitor galaxy is less ex¬ 
tended. We use a simple stellar evolution approach to compare 
our galaxy remnants to the observed dE population of the Virgo 
cluster. The aim of this study is to show the relevance of the pa¬ 
rameters, as well as their general trends with respect to transfor¬ 
mation efficiency. In section 2 we outline the numerical method 
and the different setups. Results of the different simulations are 
presented in section 3. In section 4 we compare our simulations 
with observations. In section 5 the results are discussed and in¬ 
terpreted, leading to our conclusions in section 6. 


2. Numerical methods 


There are many parameters that affect the tidal interaction of 
an infalling galaxy inside a galaxy cluster: the orbital param¬ 
eters, the structural parameters of the infalling galaxy, and the 
structural parameters of the cluster. For this study we decided 
to keep the cluster fixed and concentrate mostly on the orbital 
parameters, but including a variation of the structural parame¬ 
ters of the infalling galaxy. We used a combination of N-body 
simulations and a Monte-Carlo method to simulate the harass¬ 
ment scenario. For these simulations the direct V-body integra¬ 
tor NBO DY6 was taken (lAarsethll 19991 [2003b . with GPU accel¬ 
eration by iNitadori & Aarsethn 2012 I) and the massively p arallel 
version using MPI (ISpurzemlll999l:ISpurzem et al.ll2008l) . This 
code, called NBODY6++GPU (Wang et al., in prep.), was mod¬ 
ified to treat strong and time-varying tidal fields (see Appendix) 
and to simulate encounters with massive point-like particles for 
the harassment scenario. 

While the current version of the code can simulate up to 10 6 
particles now on a few dozen GPU nodes (order 10 5 GPU cores 
for Kepler GPUs), we used only much more moderate particle 
numbers; any unwanted two-body relaxation occurring due to 
too small a particle number has been corrected for (see Sect. 2.3). 
The code contained the SSE (single-stellar evolution) package 
of lHurlev et al] (120001) and has been combi ned with the GALEV 
spectral synthesis code dKotulla et al.ll2009l: Pang et al., in prep.). 
Although NBODY6 and its offsprings were originally designed 
for star-by-star cluster simulations, the code turned out to be very 
useful in this project, also thanks to the GPU acceleration. We 
have not resolved our galaxies star-by-star, so a single body in 
the simulation typically represents many stars. 


2.1. Cluster model 

The gravitational force of the cluster was modelled as an ana¬ 
lytical force field. An NavarroFrenkWhite profile dNavarro et al.l 
119961) was adopted to describe the density profile of the galaxy 
cluster, and th e fit paramete rs po = 3.2 • 10 -4 M o /pc -3 and 7^ = 
560 kpc from iMcLaughlinl (1 999 ) were used to reproduce the 
force field of a Virgo-like cluster with M 200 = 4.2 • 10 14 M Q 
and R 200 = 1-55 Mpc. Nbody6++GPU is using Hill’s approxi¬ 
mation to take external force-fields into account. In this approx¬ 
imation it is assumed that the galaxy is moving on an circular 
orbit around the cluster centre, which allows the Euler forces to 
be neglected. We have added some code to Nbody6++GPU to 
include the Euler forces, so that we are able to use eccentric or¬ 
bits for the galaxies around the cluster centre. An explanation of 
the equations of motions that are used is given in Appendix A. 
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The orbit of the infalling galaxy around the cluster centre was 
calculated using a leap-frog algorithm. Free parameters of the 
orbit were the apocentre and pericentre; we kept the pericentre 
fixed at 0.3 Mpc and used three different apocentres (0.5 Mpc, 
1.0 Mpc, and 1.5 Mpc). 

2.2. Encounters with other cluster members 

Flyby events of cluster members were simulated by using a 
Monte-Carlo method. A random number test was used to check 
periodically if an encounter happened or not. In the case of an 
encounter, an additional particle was injected into the simula¬ 
tion, which represents the flying-by cluster member. After that 
we waited half of the typical time between two encounter events 
and checked again whether another encounter had happened. 

The typical time between two encounters depends on the lo¬ 
cal galaxy density and the local velocity dispersion of the clus¬ 
ter. In this way the encounters with other cluster members are 
more frequent near the cluster centre. The local galaxy density 
was derived from the radial galaxy distribution of th ose galaxy 
cluste rs of the Millennium II simulation dBovlan-Kolchin et al.l 
120091) that have masses between 2.4 - 4 • 10 14 M o . We made u se 
of the corresponding semi-analytic model of lGuo et aD (1201 ll) to 
select cluster member galaxies with total (dark+stellar) masses 
of more than 0.1 M ga i, where M ga i is the initial total mass of our 
simulated infalling galaxy. The local velocity dispersion was cal¬ 
culated numerically as a solution of the spherical Jeans-equation 
for the cluster potential. 

The masses and velocities of the injected particles were cho¬ 
sen randomly. The masses following the mass function of the 
aforementioned galaxy clust ers of the Millennium II simulation 
dBovlan-Kolchin et"aLll2009l) . The minimum mass of a flying-by 
galaxy that was assumed to be relevant strong enough to disturb 
the infalling galaxy, hence is worth considering, was 0.1M ga i. 
The velocities follow a Maxwell-Boltzmann distribution, and the 
velocity orientation is chosen randomly. 

The impact parameters of the encounters were chosen by fol¬ 
lowing the geometrical cross section. Two different galaxy mod¬ 
els with different extents were used in this study. The minimum 
impact parameter for simulations with the more extended galaxy 
model was set to 8 kpc and for simulations with the smaller 
galaxy model to 4 kpc, the maximal impact parameter was set 
to 60 kpc. Owing to gravitational focusing, the distance of mini¬ 
mal approach will be closer than the impact parameter at infinity. 
To take this into account and to correct the impact parameter at 
infinity, we enlarged it by usin g the relation of the cl osest ap¬ 
proach and impact parameter of ISpurzem et al.1 ( 2009 ). A more 
detailed description on how the parameters of encounters were 
chosen is given in Appendix B. 

2.3. Galaxy models 

The galaxy models were three component models with a stellar 
disk and bulge and a dark mat t er halo created by the GalactICS 
code o f iKuijken & Dubinskil ( 1995 ). This code uses a King 
mode l dKingl 19621) for the bulge, an ex ponential Shu model dShul 
1 19691) for the disk, and an Evans model (lEvansI 19931) for the dark 
matter halo. 

Two different galaxy models were used, one with a disk scale 
length ro of 3 kpc and one with a disk scale length of 1.2kpc. 
The 3kpc model was our main model, which was used in 20 of 
24 simulations and is referred to as the ’’large model”. The 1.2 
kpc model is referred to as the ’’small model”. The large model 


Table 1. Galaxy models 


Model 

Mot 

47 s tar 

D/B 

ro 

zo 

small 

1.74 

0.12 

6.4 

1.2 

0.13 

large 

5.00 

0.35 

6.0 

3.0 

0.45 


Notes. (1) Total mass in 10 10 M©; (2) Stellar mass in 10 10 M 0 ; (3) Disk 
to bulge ratio; (4) Scale length of the disk in kpc; (5) Scale height of the 
disk in kpc 


is based on the typical va lues of dust-corrected k-b and obser¬ 
vations of Sc/Sd galaxies dGraham & Worlevll2008l) . which we 
assume to be suitable late-type progenitor galaxies. The mass-to- 
light ratios for the conversion of the k-band values were 0.46 for 
the disk and 0.82 for the bulge component. We determine these 
conversi on factors with the G ALEV evolutionary synth esis mod¬ 
els fr om Ikotulla et~akl (l2009h with Sandage-functions ( S andagel 
119861) for the star formation histories and delay times of r = 3Gyr 
for the bulge and r = 9Gyr for the disk. 

The small model is constructed to be a dwarf version of an 
Sc/Sd galaxy. The properties of both are listed in Table [TJ With 
a thickness of ro/zo = 9.2 for the small model and ro/zo = 6.7 
for the large model, both models have a comparatively thin disk 
component dBizvaev & Mitronovail20091) . The models were rep¬ 
resented by 10 5 particles, with 3 • 10 4 disk particles, 5 • 10 3 bulge 
particles, and 6.5 • 10 4 halo particles. Because Nbody6++GPU 
is a collisional N-Body code, one needs to correct the simula¬ 
tion results with respect to two-body relaxation. To determine 
the evolution under two-body relaxation, a control simulation 
was made for each of both models. In this control simulation, 
the evolution of the galaxy models without any tidal field or 
galaxy encounters was simulated, yielding the evolution of the 
model galaxy parameters th at is caused by two-body r elaxation 
alone. Following the idea of lVelazquez & White ( 19991) . one can 
correct the measured parameters from the simulations by sub¬ 
tracting the reference values of the parameter change caused by 
two-body relaxation, as determined by the control simulation. In 
that way one gets the relaxation-free change of the model param¬ 
eters. We have cross-checked this assumption, by simulating two 
of our simulations again, doubling the total number of particles, 
and it turned out that the correction works fine. The differences 
in the relaxation-free values were in the range 5% - 10%. 

2.4. Setups 

Three different orbits were used in this study: their apocentres 
were 1.5 Mpc, 1.0 Mpc, and 0.5 Mpc, and all of them have the 
same pericentre at 0.3 Mpc. The different orbits are shown in 
Fig. [TJ All simulations started at their apocentre and spanned 
a time of 5 Gyr. The time span of 5 Gyr i s similar to previ¬ 
ous s tudies (iMoore et al.lll996l Il998l 1 19991 : iMastropietro et al.l 
120051) . and it was chosen to follow the evolution of galaxies in¬ 
side a cluster environment since z ^ 0.5 . 

To probe the importance of the initial structural parameters 
of a galaxy, we ran one simulation setup with a smaller galaxy 
model. It was shown bv lYillalobos et al. ( 201 2) that the inclina¬ 
tion of the disk to the orbital plane could change the evolution of 
a galaxy in a group or cluster environment, in the way that galax¬ 
ies with an inclination could resist the tidal forces for a longer 
time. To investigate how strong this impact is, we took three dif¬ 
ferent inclinations for the galaxies on the innermost orbit: the 
disk plane parallel to the orbital plane 0°, the disk plane per- 
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Table 2. Simulation setups 


Model 

Apo 

Peri 

Inc 

— 0) 

large 

0.5 

0.3 

0 

37 

large 

1.0 

0.3 

0 

57 

large 

1.5 

0.3 

0 

73 

small 

0.5 

0.3 

0 

100 

large 

0.5 

0.3 

45 

37 

large 

0.5 

0.3 

90 

37 


Notes. (1) Used galaxy model; (2) apocentre in Mpc; (3) pericentre in 
Mpc; (4) inclination in dgre; (5) initially bound total mass in % 



[kpc] 



Time [Gyr] 


Fig. 2. Time evolution of the bound stellar mass fraction in sim¬ 
ulations with different orbital parameters of the large galaxy 
model (coloured) and with the small model (grey). The dotted 
lines represent simulations without encounters with other clus¬ 
ter members. 

All simulations started at the apocentre of their orbits. At 
this point the initial tidal radius was mostly already smaller than 
the extent of the dark matter halo and therefore only part of the 
whole mass of the galaxy was bound initially. The last column 
of Table[2] shows the initially bound mass fractions. 


Fig. 1. Orbits with a pericentre of 0.3 Mpc and an apocentre of 
0.5 Mpc (red), 1.0 Mpc (green dashed), and 1.5 Mpc (blue dot¬ 
ted). 


pendicular to the orbital plane 90°, and one orientation between 
45°. 

We made three simulations per setup, each with a different 
random number seed to get different realisations of the random 
processes. In addition, one simulation per setup was made with¬ 
out encounters with other cluster members, to study the influence 
of the tidal field alone and to investigate the stochasticity of the 
harassment scenario. Our setups are listed in Table [2] 

2.5. Bound mass 

The following algorithm was used to determine the bound par¬ 
ticles. We calculated the kinetic and potential energy of each 
particle and checked whether the absolute value of the potential 
energy was higher than its kinetic energy. In that case the particle 
was bound energetically, but it remained unclear whether it was 
bound to the galaxy or only to the cluster. Particles that were 
outside the Lagrangian points LI and L2 of the galaxy are not 
bound to the gravitational centre of the galaxy but to the cluster 
and will leave the galaxy following the attraction of the cluster. 
To distinguish between both cases, we used the tidal radius Rt , 
which is an approximation of the distance from the galaxy centre 
to the Lagrangian points LI and L2, in order to clarify whether 
a particle lies beyond the Lagrangian points or not. Because the 
tidal radius itself is only an approximation of this problem, we 
used the criteria for the discrimination that a particle has to be 
inside of 1.5 Rt around the density centre of the bulge,in order to 
be bound to the galaxy. 


2.6. Virtual photometry 


The photometric properties of different galaxy types are not ran¬ 
domly distributed in phase space, but the different galaxy types 
populate specific regions of the phase space. If the tidal transfor¬ 
mation of late-type galaxies by the harassment scenario is one of 
the main production channels of dE galaxies in galaxy clusters, 
the remnants of harassed galaxies have to match the photometric 
properties of the dE galaxies. To investigate whether the simu¬ 
lated remnants fulfil this condition, the photometric properties 
of the remnants were calculated and compared to those of dE 
galaxies in the Virgo galaxy cluster. 

To compute the brightness and the half-light radius of the 
simulated galaxies, the stellar evolution of the bulge and disk 
componen twe r e sim u lated by using the single stellar evolution 
code of [Hurley et al.l d200Qb . To do so, we adopted a Sandage- 
function (Sandage 198^) for the star formation history, following 
iGavazzi et al.l (120021) . with a delay time r = 3Gyr for the bulge 
and r = 9G yr for th e disk. As initial mass function, the Kroupa 
IMF (lKroupat200ll) in a mass range of 0.08M© - 8M© was used. 
The metallicities of the stars were calculated by using the chem¬ 
ically consistent evolution of the GALEV population synthesis 
models (iKotulla et al .1120091) . Each stellar particle was assigned 
an individual mass, age, and metallicity, and the evolution of 
the star was traced. Finally parts of the GALEV population- 
synthesis code (IKotulla et al.l2009l) were used to convert the stel¬ 
lar properties i nto the r-band brightness o f the stars by using the 
BaSeL-library dLeieune et al .11 19971 1 19981) . 

Furthermore, it was assumed that the star formation ended 
abruptly when the galaxy entered the cluster. While it is ex¬ 
pected that a low-mass galaxy loses most of its gas relatively 
quickly by ram pressure ( Gunn & Gottll 19721; iBoselli et al.ll2008l) 
- eve n outside the virial radius when entering along filaments 
dBahe et al.l [201 3?) - this is certainly a simplified approach. 
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Models employing an instantaneous removal of the extended gas 
reservoir of galaxies that become satel lites led to an overabun - 
dance of passive galaxies in clusters dWeinmann et al.l [20 id) : 
a recent study combining observations and simulations instead 
points to a pictu re in which ’’satellite quenching is delay ed-then- 
rapid” (I Wetzel et al.1 (2013h . On the other hand, the relevance of 
pre-processing (IVillalobos et al.ll2012l) may imply an earlier de¬ 
cline of the star formation rate. Since we are primarily interested 
in simulating the effects on the stellar structure, we consider the 
star formation rate to be zero from the starting point of our sim¬ 
ulation, 5 Gyr ago. 

To compute the total brightness and the half-light radius of 
the galaxies, projections of the galaxies were created under dif¬ 
ferent projection angles 0 e [0°, 22.5°, 45°, 67.5°, 90°] along the 
x and y-axes. We extracted radial pr ofiles by fitting ellipses t o 
the isophotes following the formula of lBertin & Arnouts ( 199 6). 


As total brightness of the galaxy, the brightness inside of an 
elliptical aperture of two Petrosian semi-major axes was used 
(|Petrosianlll976l: iLisker et ai1l2007l) . The major axis of the el¬ 
liptical profile which encloses half of the total luminosity was 
adopted as the half-light radius of the galaxy. From the total 
brightness and the half-light radius one can compute the effec¬ 
tive surface brightness. 

For the comparison with o bservations, dat a fro m the Sloan 
Digital Sky Survey (SDSS) (lAdelman-McCarthv et al. 2006 ) 
were used. The sample consists of 864 members and possible 
members of the Virgo cluster, classified as dE, dEdi, SO, Sc, 
Sd, Sm, and Im galaxies based on the Virgo Cluster Catalogue 
(iBinggeli et al.l ll98~5h down to the completeness limit of a B- 
band apparent magnitude 18 and with velocities below 3500 
km/s. The measurements of r-band brightn ess, surface bright¬ 
ness, and half-light radius w ere t aken from ILisker et al.l (120071) . 
iJanz & Liskerl (120081120091) . and lMever et al.l (120141) . The sam¬ 
ple comprises 449 dE , 41 dEdi, 73 SO, 114 Sc/Sd, 26 Sm, and 
161 Im galaxies. 


3. Results 

In the following section, we want to analyse the impact of or¬ 
bital parameters, the structural parameters of a galaxy, and the 
inclination of a galaxy. In Table [3l the initial properties of the 
galaxy models, as well as the properties of the galaxy remnants 
of the different simulations, are listed. The values of Tableware 
already relaxation-corrected as explained in section 2.3. 

3.1. Effect of orbital parameters and encounters 

As expected, the orbital parameters were essential for the evolu¬ 
tion of a galaxy inside of a galaxy cluster environment. The evo¬ 
lution of the bound stellar mass fractions in the simulations with 
different apocentres (0.5Mpc, l.OMpc, and 1.5Mpc) is shown in 
fig. [2 On the farthest orbit, the stellar mass loss was between 4% 
and 6%, whereas on the innermost orbit, the stellar mass loss was 
24% to 50%. On the intermediate orbit, it reached 8% to 13%. 
The scatter of the mass loss in the different setups is a result of 
the stochastical nature of the harassment scenario. 

Figures [3] and|4] show the evolution of the galaxies in the first 
and second runs of the setup with the closest orbit by a series of 
snapshots with intervals of 625 Myr. In the first run of this setup, 
no encounter with another cluster member resulted in additional 
mass loss. As opposed to that, in the second run at a very early 
stage of the simulation, an encounter pulled nearly 5% of all stel¬ 
lar particles out of the disk and disturbed the system strongly. 


This resulted in a different mass loss history and different shape 
of the remnant. For a better understanding of the mass loss his¬ 
tories and evolution of the galaxies, it is reasonable to separate 
the effects of the global tidal field and close encounters, as we 
do in the following. 

3.1.1. Mass loss on different orbits without encounters 

Using the simulations without encounters we could follow the 
evolution of the infalling galaxy introduced only by the global 
tidal field of the cluster. This evolution is a kind of minimal ef¬ 
fect that each galaxy experiences if it falls on a given orbit into 
a cluster. That effect could be amplified by encounters, but not 
weakened. 

The stellar mass loss by the global tidal field (dotted lines 
in Fig. [2]) started when the galaxy passed the pericentre the first 
time. At the pericentre the tidal forces reached their maximum 
strength and therefore the tidal radius reached a local minimum. 
The local minimum of the tidal radius is equivalent to a local 
minimum in the bound mass fraction. The small tidal radius at 
the pericentre enabled the growth of tidal tails. Along these tidal 
tails, stars also left the galaxy more than one Gyr after the galaxy 
passed the pericentre. The smaller the apocentre of an orbit, the 
more frequently the galaxy passes the pericentre and the smaller 
the tidal radius averaged over the simulation time. The tidal force 
becomes more efficient in removing material the smaller the av¬ 
eraged tidal radius gets relative to the averaged half-mass radius 
of the stellar component of the galaxy. Figure \5\ shows that for 
the galaxies without an inclination, the ratio of the averaged tidal 
radius along the orbit and the averaged half-mass radius of the 
stellar component is an adequate way to quantify the efficiency 
of tidal forces. The mass loss caused by the persistent global 
tidal forces of the cluster’s background mass represents the min¬ 
imum mass loss by the harassment scenario. 

A huge part of the dark matter halo was also stripped at the 
first pericentre passage. For the large model, the loss of initially 
bound total mass was close to or more than 50%, for the small 
model this loss was in most cases only 6%. 

3.1.2. Mass loss on different orbits with encounters 

Within one setup, the runs show a large scatter in their stellar 
mass loss. This scatter results from different encounters with 
cluster members. Not every encounter has the same impact, it 
depends on the minimal distance, the mass ratio and the relative 
velocity of both galaxies. To describe the impact of encounters 
in a similar way to how we described the impact of the global 
tidal field before, we calculated the tidal radius of the galaxy 
with respect to the flying-by objects at the moment of closest 
approach. This tidal radius is expressed in units of the half-mass 
radius of the stellar component of the galaxy model. This quan¬ 
tity shows whether a flying-by cluster member could have any 
influence on the galaxy. However, one also needs to take the rel¬ 
ative velocity of both galaxies into account, because the slower 
a galaxy is as it fies by, the higher the energy and angular mo¬ 
mentum transfer from this object. Figure[6| shows the distribution 
of the above properties for all flying-by galaxies in our simula¬ 
tions. The symbol colour indicates the approximate energy that 
the flying-by cluster members inject into the stellar component 
of the galaxy, in units of the total energy of the stellar compo¬ 
nent. To calculate this deposited energy, we used the approxi¬ 
mation of ISpitzei] (19581) for impulsive high-speed encounters. 
This is a simplification of the real galaxy-galaxy encounter, but 
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Fig. 3. Snapshots from the first simulation run of the setup with the closest orbit, i.e. an apocentre of 0.5 Mpc and pericentre of 0.3 
Mpc. The timestep between snapshots is 625 Myr; they start at the upper left at 0 Gyr and reach the lower right after 5 Gyr. Each 
snapshot image shows the galaxy face-on and from two perpendicular edge-on viewpoints. In the upper right corner of each image, 
the actual position on the orbit around the cluster centre is marked. The red dots are bound disk particles, the green dots are bound 
bulge particles, and the black dots are unbound stellar particles. 


it is accurate enough to give an impression of the strength of 
the encounter. The figure illustrates that only a small number of 
all encounters with other cluster members fulfilled both crite¬ 
ria: to have a strong influence and be slow in relative velocity. 
Consequently, only a few encounters result in significant energy 
injection and additional mass loss. Therefore several encounters 
with cluster members are needed to have a strong impact on the 
infalling galaxy. This small fraction of particularly strong en¬ 
counters with cluster members is also why the mass loss curves 


of several simulations are close to the mass loss curves of sim¬ 
ulations without encounters. The different evolution of a galaxy 
with and without strong encounter is illustrated in Figs. 3 and 4 
for the same simulation snapshots. The total number of encoun¬ 
ters and the number of strong encounters (with A E > 10 -2 ) in 
each simulation run is given in columns (6) and (7) of Table [3] 
The total number of encounters in simulations with the small- 
galaxy model is two times higher than for simulations with the 
large model on the same orbit, because for the encounter proba- 
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Fig. 4. Like Fig. El but for the second simulation run of the setup with the closest orbit. 


bility, we take cluster members with masses above 0.1M ga i into 
account(M ga i is the initial mass of the galaxy model), therefore 
more galaxy encounters occur to the smaller and lighter model. 
But it is also harder to affect the small model by an encounter. 
It has to come closer to the galaxy to strip the disk stars of the 
smaller model away, and the probability of such close encounters 
declines quadratically with the impact parameter. 

The scatter of mass loss within a given setup declines with 
the increasing apocentre of the orbit. The reason is that, in setups 
with larger apocentre, the galaxies have spent a smaller amount 
of time in regions of the cluster with a high galaxy density, and 
therefore fewer encounters occurred. 


3.1.3. Transformation of shapes and kinematics 

The top panel of Fig. [7] shows the thickness of the galaxy rem¬ 
nant against the ratio of the average tidal radius to the aver¬ 
age half-mass radius. The thickness of the galaxies is measured 
as the ratio of the square root of the quadratic mean height of 
stars and the square root of the quadratic mean radius of stars 
{c/a = a/(z 2 )/ >/(r 2 )). When averaging over the runs of each 
setup, we find that the stronger the tidal force, the stronger the 
average thickening of the galaxies. At a fixed ratio of averaged 
tidal radius to averaged half-mass radius, the thickening was 
increased by the individual encounters with cluster members. 
Within a given setup, the stronger the stellar mass loss induced 
by such an encounter, the stronger the induced thickening. But 
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Table 3. Simulation results 


Setup 

Galaxy properties 

Apo. 

Inc. 

Model 

Enc. 


#enc #enc 

'aE>10~ 2 

Star 

Tot 

Rl/2 

< z > 

c/a 

^max 

cr 

^max / U" 

Properties of the initial models 

- 

- 

large 

- 

- 

- 

- 

100 

100 

5.00 

0.37 

0.08 

65 

25 

2.60 

- 

- 

small 

- 

- 

- 

- 

100 

100 

2.00 

0.16 

0.09 

89 

39 

2.28 

Properties of the galaxy remnants 

0.5 

0° 

large 

y 

3.7 

15 

0 

77 

46 

3.40 

0.60 

0.16 

48 

39 

1.22 

0.5 

0° 

large 

y 

3.7 

20 

1 

66 

32 

3.10 

1.05 

0.32 

30 

36 

0.84 

0.5 

0° 

large 

y 

3.8 

18 

1 

55 

30 

2.50 

1.12 

0.43 

29 

43 

0.68 

0.5 

0° 

large 

n 

3.7 

0 

0 

75 

46 

3.50 

0.59 

0.17 

48 

37 

1.28 

1.0 

0° 

large 

y 

5.3 

14 

0 

92 

39 

3.90 

0.77 

0.15 

46 

41 

1.12 

1.0 

0° 

large 

y 

5.5 

14 

0 

89 

37 

3.60 

0.85 

0.19 

42 

39 

1.07 

1.0 

0° 

large 

y 

5.3 

13 

1 

87 

35 

3.90 

1.41 

0.21 

40 

36 

1.12 

1.0 

0° 

large 

n 

5.3 

0 

0 

92 

40 

3.90 

0.67 

0.11 

47 

39 

1.20 

1.5 

0° 

large 

y 

7.8 

11 

0 

96 

40 

4.20 

0.59 

0.10 

52 

36 

1.45 

1.5 

0° 

large 

y 

8.0 

9 

0 

94 

37 

3.90 

0.61 

0.10 

52 

37 

1.39 

1.5 

0° 

large 

y 

7.8 

5 

0 

96 

38 

4.10 

0.58 

0.10 

53 

36 

1.48 

1.5 

0° 

large 

n 

7.8 

0 

0 

96 

40 

4.10 

0.53 

0.08 

53 

32 

1.64 

0.5 

0° 

small 

y 

8.4 

42 

0 

100 

94 

2.00 

0.17 

0.09 

87 

41 

2.00 

0.5 

0° 

small 

y 

8.4 

44 

0 

98 

77 

2.00 

0.43 

0.19 

81 

41 

1.99 

0.5 

0° 

small 

y 

8.4 

45 

0 

100 

94 

2.00 

0.17 

0.09 

86 

41 

2.12 

0.5 

0° 

small 

n 

8.4 

0 

0 

100 

94 

2.00 

0.16 

0.08 

87 

41 

2.14 

0.5 

45° 

large 

y 

3.6 

13 

0 

80 

46 

3.60 

1.11 

0.41 

50 

37 

1.33 

0.5 

45° 

large 

y 

3.4 

18 

1 

28 

14 

2.80 

1.37 

0.49 

22 

29 

0.76 

0.5 

45° 

large 

y 

3.8 

21 

0 

72 

41 

3.20 

1.43 

0.54 

43 

36 

1.20 

0.5 

45° 

large 

n 

3.7 

0 

0 

80 

46 

3.50 

1.11 

0.42 

50 

37 

1.33 

0.5 

90° 

large 

y 

3.4 

19 

0 

88 

51 

4.00 

0.81 

0.29 

59 

31 

1.93 

0.5 

90° 

large 

y 

3.4 

21 

1 

77 

43 

3.60 

1.47 

0.34 

42 

36 

1.17 

0.5 

90° 

large 

y 

3.4 

20 

2 

62 

32 

3.10 

1.86 

0.55 

34 

37 

0.91 

0.5 

90° 

large 

n 

3.4 

0 

0 

89 

49 

4.00 

0.89 

0.33 

60 

31 

1.96 


Notes. (1) apocentre of the orbit in Mpc, the pericentre was always fixed to 0.3Mpc; (2) inclination of the orbital plane and the disk; (3) used 
galaxy model; (4) enc. indicates if a simulation is using encounters with cluster members or not; (5) ratio of the averaged tidal radius and the 
half-mass radius; (6) and (7) are the number of galaxy-galaxy encounters in this simulation and the number of galaxy-galaxy encounters with an 
approximated energy transfer greater than 10 -2 ; (8) and (9) fraction of the initially bound stellar and total mass; (10) half-mass radius of the galaxy 
in kpc; (11) mean height of stars above the disk in kpc; (12) axis ratio of the stellar component; (13) maximum rotation velocity in km/s; (14) 
velocity dispersion in km/s; (15) ratio of peak velocity and mean velocity dispersion. All the values have already been corrected for relaxation as 
explained in section 2.3 
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Fig. 5. Bound stellar mass fraction versus the averaged tidal 
radius in units of the averaged half-mass radius of the model 
galaxy. Open symbols represent simulations without encounters 
with other cluster members. 


also without additional mass loss, encounters were able to dis¬ 
turb the disk and induce a small thickening of the galaxy. The 


Fig. 6 . Relative velocity of the flying-by cluster members versus 
the tidal radius with respect to the flying-by cluster members in 
units of the half-mass radius of the model galaxy at the moment 
of closest approach. The colour scale denotes the approximate 
energy that the flying-by cluster members inject into the stellar 
component of the galaxy, in units of the total energy of the stellar 
component. 
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Fig. 7. Top: Axis ratio of the stellar component of the galaxy 
remnants versus the averaged tidal radius in units of the aver¬ 
aged half-mass radius. Bottom: Ratio of the rotational peak ve¬ 
locity and the mean velocity dispersion of the galaxy remnants 
versus the averaged tidal radius in units of the averaged half¬ 
mass radius of the stellar component. Both panels only contain 
simulations without any inclination between orbital plane and 
galactic disk. Open symbols represent simulations without en¬ 
counters with other cluster members. 


bottom panel of Fig. [7] shows the ratio of maximal rotation ve¬ 
locity and the mean velocity dispersion of the stellar component. 
The galaxy model started with a rotation-dominated stellar com¬ 
ponent, but the rotational support decreased with time by mass 
stripping and disturbances. Even if only small or no additional 
stellar mass loss happened - as in the case of the simulations 
with an apocentre of 1.5 Mpc - the ratio of peak velocity and 
mean dispersion could be decreased by perturbations by ~ 10% 
or more. In the two simulations where the stellar mass loss was 
above 30%, the ratio of the peak velocity and the mean velocity 
dispersion drops below one and the systems become dominated 
by pressure. 

3.2. Effect of initial structural parameters 

Changing the initial extent and the mass of a galaxy model is an¬ 
other way to change the ratio of the averaged tidal radius and the 
averaged half-mass radius of the stellar component while keep¬ 
ing the orbital parameters fixed. For a better comparison we plot 



Time [Gyr] 

Fig. 8. Bound stellar mass fraction in the simulations with differ¬ 
ent orientations of the galactic disk. The dotted lines represent 
simulations without encounters with other cluster members. 

the results of the simulations with the small model, together with 
those of section 3.1, where we used the large model in figures O 
\5\ and El It illustrates how important the initial structural param¬ 
eters of a galaxy model are. They are as relevant as the orbital 
parameters. 

In the simulation without encounters, the small model lost no 
stellar mass, experienced no thickening, its ratio of peak veloc¬ 
ity and mean dispersion decreased only by ~ 6% and only ~ 6% 
of its dark matter mass got lost. The small model also had a 
smaller cross section for encounters with other cluster members, 
which made it hard to be influenced by such encounters. Only 
one encounter in one of the runs was strong enough to disturb the 
galaxy, and it took only 2% of the stellar mass away. However, 
this single encounter increased the thickness of the disk by a 
factor of two and decreased the ratio of maximal rotation veloc¬ 
ity and mean dispersion by an additional 8%. These results fit 
well with the finding of sections 3.1 and 3.2 that the extent of a 
galaxy has to be close to its tidal radius and the galaxy has to be 
in a region with high galaxy density to get tidally transformed 
efficiently. 

3.3. Effect of inclination 

The inclination of the galactic disk relative to the orbital plane 
is an additional parameter that influences the efficiency of the 
tidal transformation of galaxies. We ran simulations for galax¬ 
ies with different inclinations (0°, 45°, and 90°) on those of 
our orbits that are farthest inside. An inclination is expected 
to de crease the efficiency o f mass stripping by the global tidal 
field (I Villalobos et al.|[2Q12h . but in general it is not expected to 
change the effect of close encounters. Figure [8] shows the evo¬ 
lution of the bound stellar mass fraction in the simulations with 
different inclinations and the closest orbit. As expected the mass 
loss induced by the global tidal field alone (dotted lines in Fig. 
[8]) decreases with increasing inclination. In the case of a perpen¬ 
dicular orientation of the disk (90°), the stellar mass loss induced 
by the global tidal field was reduced by half compared to the case 
where the disk is parallel to the orbital plane (0°). The additional 
mass loss caused by encounters with other cluster members still 
shows a large scatter. 

Figure [9] shows the thickness and the ratio of peak velocity 
to the mean velocity dispersion against the bound stellar mass 
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Fig. 9. As in Fig. [7] but for simulations with an inclination be¬ 
tween orbital plane and galactic disk. 

fraction. Galaxies with the same stellar mass loss show stronger 
thickening if they do not move parallel to their disk. Under the 
influence of the global tidal field and without encounters, the 
galaxies with an inclination thicken more. Especially with an in¬ 
clination of 45°, the galaxy thickens by a factor of five, whereas 
in the case of 0°, the galaxy thickens only by a factor of two. 
The simulations without encounters show that the galaxy with 
a perpendicular orientation of the disk to the orbital plane could 
preserve its rotational support much better than the galaxies with 
other inclinations. Nevertheless, under the influence of encoun¬ 
ters, a galaxy with perpendicular orientation of the disk also be¬ 
came pressure-dominated. This highlights the stochastical nature 
of the harassment scenario again. 

One has to stress that the effects an inclination has on the 
interaction with the global tidal field of a cluster are diverse: it 
could reduce the mass loss, amplify the thickening, and reduce 
the kinematical transformation. However, it seems that the in¬ 
fluence on each of these quantities differs between the different 
orientations. Therefore, to describe the evolution of galaxies in 
a cluster environment, one has to pay attention to the different 
possible orientations of the disk. 

4. Comparison with observations 

Figure [T0| shows the distribution of different galaxy types of 
the Virgo cluster dLisker et al.ll2QQ7b iJanz & Linke r 2008, 20091: 


iMever et ai1l2014l) and the distribution of our initial galaxy mod¬ 
els (green), as well as their remnants (red) in the (yu) eff - r e ff (top 
panel) and in the M r - (yu) eff (bottom panel) plane. The points 
for the initial and remnant galaxies are the mean values averaged 
over different projection angles (edge-on = 90°, 67.5°, 45°,22.5°, 
0° = face-on) along the x-axis and along the y-axis; the error 
bars span the range from the minimum to the maximum value of 
these projections of the galaxies. The small scatter of the initial 
models originates in the different random-number seeds that are 
used to assign masses and ages to the stellar particles. The large 
error bars in the surface brightness of the initial models stems 
from the fact that the surface brightness of the edge-on projec¬ 
tion of the model is nearly 2 mag/arcsec 2 brighter than the face- 
on projection. The point clouds of initial and remnant galaxies 
that belong to simulations with the large model can be identified 
by their clear separation in effective radius and magnitudes in 
comparison to the point clouds of the smaller model. 


Compared to the distri bution of the Sc/Sd galaxies in the 
Virgo cluste r as given by Lisker e t al.1 (l2QQ7h . Uanz & Liskerl 
(120081 120091) . and Meyer et al.l (120141) . the i nitial large models 
that ar e based on the Sc/Sd parameters of iGraham & Worlevl 
( 2008 ) tend to have larger radii. The initial small models were 
smaller than the typical extent of the Virgo Sc/Sd galaxies, which 
have a mean effective radius of = 2.63 kpc. The point cloud 
of remnants of the large model overlaps partially with the more 
extended dEs and dEdis. But on average the remnants of the 
large model tend to have radii that are too large to fit with the 
distribution of dEs well. Their magnitudes are close to -18 mag 
in r, which is often used as the threshold for classifying bright 
dwarf galaxies, but in general they are too bright for ordinary dEs 
or dEdis. This trend indicates that the stripping of stellar mate¬ 
rial, hence the truncation of the disk and the dimming of light 
in simulations with the large model, was not efficient enough to 
produce typical dE or dEdi galaxies. However, the direction of 
evolution of the large model implies that harassment of galax¬ 
ies could potentially produce dE and dEdis if the effects were 
amplified in some way. 


The remnants of the small model are all placed inside the 
distribution of bright dE and dEdi galaxies. But the dimming of 
light of these remnants was only caused by the truncation of the 
star formation at the beginning of the simulation, not by strip¬ 
ping of stellar matter. We have to stress that it has been shown 
before that these galaxies were nearly unaffected by the envi¬ 
ronment, so these remnants are still disk galaxies with a small 
thickness (c/a = 0.09 - 0.19) and are strongly dominated by ro¬ 
tation (v/cr = 1.99 - 2.14). They could not be used to explain the 
distribution of dE galaxies because they did not undergo a mor¬ 
phology transformation. Nonetheless, we can imagine that if the 
effect of the tidal interaction was amplified, then the direction of 
evolution of these galaxies would be similar to that of the large 
model, and they would undergo a morphological transformation 
and end up as dE galaxies. 


This comparison shows that there is already a partial overlap 
of harassed galaxies with dE/dEdi galaxies after an evolution of 
5 Gyr in a cluster environment. From these results one can ex¬ 
pect that if the environmental effects were amplified, the harass¬ 
ment scenario could be one production channel of dE and dEdi 
galaxies. 
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Fig. 10. Top: Distribution of Virgo galaxies, as well as the initial galaxy models and remnants in the (yu) eff - r e ff plane. Bottom: 
Distribution of Virgo galaxies, as well as the initial galaxy models and remnants in the M r - (yu) eff plane. The points for the initial 
and remnant galaxies are the mean values averaged over different projection angles (edge on = 90°,67.5°,45°,22.5°, face on = 0°) 
along the x-axis and along the y-axis. The error bars span over the range from the minimum to the maximum value of the projections. 
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5. Discussion 

5.1. Orbital parameters 

By varying the apocentre of the orbit of an infalling galaxy be¬ 
tween 0.5 Mpc and 1.5 Mpc at a fixed pericentre of 0.3 Mpc, 
it turned out that the tidally induced morphological transforma¬ 
tion of low-mass late-type galaxies to bright dE galaxies became 
more efficient, the closer the orbit of the galaxies was. With de¬ 
creasing apocentre, the ma ss loss, the thickness, and the pressure 
s upp ort increased (also see lMastropietro et ai]l2005l: [Smith et al.l 
120101) . However, the comparison of the galaxy remnants with 
observations of dE galaxies shows that even on the innermost 
orbit, the environmental effects were still not efficient enough to 
transform the infalling galaxies into typical dE galaxies within 
an evolution time of 5 Gyr. In general, the remnants are still too 
extended and too bright, which indicates that a stronger environ¬ 
ment would be necessary to strip more material and complete 
the transformation, or the progenitor galaxies would have to be 
late-type disks of lower luminosity and size than our model (see 
[Graham et al.lll2003 ). 

If one takes dynamical friction into account, galaxies that en¬ 
ter a cluster will sink deeper into the cluster core. We neglected 
this effect in our simulations and probed what happened to galax¬ 
ies that did not end up near the cluster centre. This is a relevant 
question because one observes dE galaxies not only in cluster 
cores but also in the outskirts of clusters, so a main production 
process of dE galaxies would also have to work in these regions 
and not just in the very cluster centre. 

On our more extended orbit (apocentre 1.5 Mpc) in the out¬ 
skirts of the cluster, the tidal forces of the cluster were not strong 
enough and the galaxy encounters not frequent enough to ha¬ 
rass the galaxies. On this orbit, the galaxies passed the pericentre 
only once in 5 Gyr, and the mass loss, as well as the thickening 
of th e galaxies , was marginal. These results agree with the study 
of [Smith et al.l ( 2010 ). who found in their simulations that low- 
mass dwarf irregular galaxies that enter a galaxy cluster from 
its present-day virial radius were mostly unaffected by the tidal 
forces of the cluster or by interactions with outer cluster mem¬ 
bers. Only on those orbits that are very close to the cluster centre 
did I Smith et al.l find the tidal interactions to be efficient to trans¬ 
form dwarf irregular galaxies. They conclude that galaxies that 
are transformed in a harassment scenario have to have entered a 
galaxy cluster a long time ago to reach the innermost regions of 
a galaxy cluster and get transformed. 

iDe Lucia et al.l (120121) have shown by us ing the data of the 
Millennium simulation dSpringel et al.ll2005l) that galaxies with 
halo masses below 10 n M o that are in a region closer than 0.5/?^ 
today became satellites of bigger systems 9 Gyr ago (median 
value). Nine billion years are nearly double the time that we sim¬ 
ulated. However, we could expect from our simulation results 
that a galaxy that spends such a long time in an environment like 
a cluster core could complete the transformation to a dE galaxy. 

One point that needs to be taken into account when simulat¬ 
ing such a long time period is that the cluster itself will evolve on 
this time scale. A galaxy that enters a cluster nine billion years 
ago enters a protocluste r enviro nm ent w ith a mass that is more 
like present-day groups. [Villalobos et al l (2 012 ) have done sim¬ 
ulations of the evolution of galaxies in a group-like tidal field to 
study the relevance of preprocessing of galaxies in groups. They 
find that the group tidal field can induce a morphological trans¬ 
formation (also see lMaver et al .1120011) . Taking the preprocessing 
in galaxy groups into account could thus also be an answer for 
how harassed galaxies could populate the outskirts of clusters, 


because not all protocluster members will end up in the cluster 
core. 

iLisker et al i (2007 ) mention that the different infall times and 
the associated cluster centric distances could be an explanation 
of the origin of different dE subpopulations like dEdis and could 
give a natural answer to the question of why the different sub¬ 
populations of dEs follow a morphology-density relation. In this 
sense one could interpret dEdis not only as galaxies in an un¬ 
completed transformation phase (iMastropietro et al.l 120051) but 
also as a separate subpopulation that evolved largely in a dif¬ 
ferent environment than ordinary dEs. In the same sense, all the 
different subpopulations of dEs may have evolved in parallel in 
different environments and would then be the present-day prod¬ 
uct of this parallel evolution (iLisker et al.ll201 3). 

The next step in clarifying the relevance of the harassment 
scenario as a production channel of dE galaxies would have to 
be studies that simulate a longer period in time and take the evo¬ 
lution of the environment itself into account. 

5.2. Structural parameters 

We have probed the influence of structural parameters of the pro¬ 
genitor galaxy by using two different models: one with a small 
and one with a large disk component, which differ by a factor 
2.5 in effective radius. While the size of the large model be¬ 
comes smaller by 26% on average, the size of the small model 
remains unchanged (see Table [3]). The small model was sitting 
so deep inside its tidal radius that it remained nearly untouched 
by the global tidal field. The slight increase in the half-light ra¬ 
dius (Fig. is due to the older stellar population of the bulge 
component. 

Small galaxies are much harder to affect by the tidal force 
of the cluster or by encounters with other cluster members, in 
agreement with the findings of [Moore et al. ( 1999 ). They stud¬ 
ied the efficiency of the harassment scenario for low and high 
surface brightness galaxies and found that high surface bright¬ 
ness galaxies are very resistant to environmental effects. 

A galaxy like the small model has to be on an orbit closer 
to the cluster centre to undergo a similar evolution as the large 
model. To reach the same ratio of mean tidal radius and half¬ 
mass radius, the small model would have to be about two times 
closer to the cluster centre than the closest orbit we used for the 
large model if we assume that the enclosed cluster mass stays 
the same. If one takes the radial mass profile of the cluster into 
account, one can in fact estimate that the small model would 
have to be about seven times closer to the cluster centre. This 
means it would have to be on an orbit with an apocentre closer 
than 100 kpc to the cluster centre. This emphasizes again that 
the initial structural parameters of the progenitor galaxy are as 
relevant as the orbit of the galaxy. 

5.3. Inclination 

The inclination of the disk of an infalling galaxy has been found 
to be an important parameter for the influence of the cluster’s 
tidal field. The induced stellar mass loss was reduced by half 
when the disk was perpendicular to the orbital plane as com¬ 
pared to a parallel orientation on the same orbit. On the other 
hand, the thickening of galaxies with an inclination was stronger 
than for parallel oriented galaxies. Especially with an inclina¬ 
tion of 45°, the thickening induced by tidal forces was more 
than twice as large as the parallel oriented galaxy. Additionally, 
it turned out that galaxies with a perpendicular orientation could 
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shield themselves better against kinematical transformations by 
the ti dakfie l d. 

[Villalobos et al.1 (2012 ) have found in their simulations of 
the evolution of galaxies in a group-like tidal field that galax¬ 
ies that are initially face-on or retrograde-oriented lost less mass 
and preserved their disk structure and kinematic for a longer 
time than prograde galaxies. Q They conclude that this could be 
an explanation for disk gal axies in the inner regio ns of groups 
and clusters. We agree with I Villalobos et all (120121) that galaxies 
with an inclination experienced less mass loss and kinematical 
transformations. 

We speculate that the effect of inclination could be an origin 
for dEdi gala xies. dEdis appear to be less round than ordinary 
dEs dLisker et al.ll2006bb . but with a thickness of 0.35 - 0.4, they 
are somewhat thicker than typical late-ty pe di sk g alaxies at sim¬ 
ilar stellar mass (cf. ISanchez-Janssen et al.ll2010l) . The rotation 
parameter v max /cr of dEdi galaxies is smaller than for disk galax¬ 
ies bu t higher than for ordinary dE galaxies dToloba et al.ll201 ll 
120141) . We found that our simulated galaxies with an inclination 
of 45° or 90° thicken more and reach the regime of dEdi galaxies. 
Galaxies with a disk perpendicular to their orbital plane (90°) 
could shield their kinematics much better against the influence 
of the tidal field. Both observations together lead to the specu¬ 
lation that disk galaxies that enter a cluster at a high inclination 
45° - 90° have a good chance of ending up with a thick disk with 
high v max /cr values. In that case, harassed high-inclination disk 
galaxies may be an origin for dEdi galaxies. 

6. Conclusions 

Harassment of galaxies is a complex process. Its efficiency de¬ 
pends on a combination of many parameters, such as the orbit, 
the inclination, and the structural parameters of the progenitor 
galaxy, as well as on the properties of the galaxy cluster itself. It 
is likely that harassment is one of the most important engines of 
galaxy evolution inside galaxy clusters, and is often referred to 
in past studies. However, this process would have had to start at 
a very early epoch in protocluster environments to complete the 
transformation by today. To fully understand the overall impor¬ 
tance of harassment, more simulations that reproduce the com¬ 
plicated interaction between galaxies and their environment are 
needed. 


Appendix A: Equations of motions 

The reference frame that was used to describe the equations of 
motion of a particle in our simulation was centred on the density 
peak of the bulge. The galaxies themselves moved on an eccen¬ 
tric orbit around the cluster centre. Additionally, we rotated our 
coordinate system in such a way that the x-axis of the reference 
frame always pointed away from the centre of the cluster. This 
means that our coordinate system is not an inertial system, but 
an accelerated and rotated one. 

The equation of motion of a particle is 

? = /gal + /cl + «0 + fee + fco + fe u- (A.l) 

The first two terms are the gravitational forces of the galaxy / ga i 
and the cluster / c j, which acts on a particle. The gravitational 
force of the galaxy was calculated by direct summation of the 
two-body interaction of the particles, and the gravitational force 
of the cluster was calculated by an analytical force field as de¬ 
scribed in section 2.1: 


/gal = -VOg a i (A.2) 

/d = -VO d . (A.3) 

The third term of equation A. 1 is the force of inertia «o- It has the 
same strength but inverse direction as the force that accelerates 
the galaxy on its orbit around the cluster centre. The distance of 
the cluster centre and the origin of ordinates is ro. The direction 
of the force of inertia points away from the cluster centre, which 
means that in our rotating reference frame along the direction of 
the x-axis, that 


«o 


d<b d 


dr c i 


(A.4) 


Since it is a rotating reference frame, we have to take the cen¬ 
trifugal / ce , Coriolis / co and Euler force / eu into account. The 
orbital plane of the galaxy within the cluster is chosen so that the 
angular velocity of the motion of the galaxy around the cluster 
centre is pointing in the z direction. This angular velocity is the 
same as the angular velocity that is used to rotate the reference 
frame: 
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1 Our simulations only included prograde rotation. 


/ce =-(o x ((o x r) = co 2 [xe x + yey) (A.5) 

/co = -2cj x r = 2co (ye x - xe (A.6) 
/eu = ~(oxr = (b (ye x - xe y ) . (A.7) 

After an expansion of the cluster potential in a Taylor series to 
the first order around the origin of ordinates, the equations of 
motions are 


®gal + I 


dx 


+ I" - 


3 2 Or 


dr c l 2 
dOd 


x + 2 coy + (by 


cl¬ 

oy 


< 9 ® d 


dr d 


dr c l 

z 

ro 


— I y — 2 cox - (bx 
r o 


(A.8) 
(A.9) 
(A. 10) 


If one assumes a circular orbit for the galaxy tb = 0, these equa¬ 
tions of motion simplify to the Hill’s approximation. 
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Appendix B: Properties of encounter galaxies 

The typical time between two encounters depends on the lo¬ 
cal galaxy density n(r) and the local velocity dispersion <x(r) of 
the cluster, which are functions of the cluster centric distance r. 
We take all those encounters into account that pass the infalling 
galaxy within a distance of p max = 60kpc : 


T = (2 n(r) cr(r)) -1 


(B.l) 


where 2 is the geometrical cross section of such a galaxy-galaxy 
encounter E = 717 ?^. . The local galaxy density was derived by 
assuming a galaxy distribution that follows a power law: 


dN_ 

dr 


N 0 A 


(B.2) 


The parameters No = 154 and = 0.6 were determined by a fit 
to the radial galaxy distri butio n of t h ose g alaxy clust ers of the 
Millennium II simulation dBovlan-Kolchin et al .1120091) that have 
a mass between 2.4 - 4 • 10 14 M 0 . 

Cluster members with total masses down to 0.1 M ga i were 
taken into account, where M ga i was the mass of the infalling 
galaxy. The total masses of the galaxies were taken from the 
corresponding semi-analytic model of dGuo et alJl201lh . The lo¬ 
cal velocity dispersion was calculated numerically as a solution 
of the spherical Jeans-equation for the cluster potential. The oc¬ 
currence of an encounter was simulated by randomly placing the 
particle that represents the perturber galaxy on a sphere with a 
radius of 200kpc around the infalling galaxy. The velocities of 
the flying-by galaxies followed a Maxwell-Boltzmann distribu¬ 
tion. The velocity vector of the flying-by galaxy pointed to the 
centre of the sphere, but offset by the impact parameter of the 
encounter. The impact parameter was chosen randomly based 
on the geometrical cross section, so that the probability of an 
impact parameter p e \p m i n , Pmax ] is given by 


m = 2 



(B.3) 


One also has to consider that the flying-by galaxy will by 
focused on the infalling galaxy by gravitation. We corrected the 
offset to take this g ravita tional focusing into account by follow¬ 
ing ^purze^^^Jd^^l) . 

The masses of the perturber galaxies were chosen randomly 
by following the mass distribution of the a forem ention ed galaxy 
cluste rs of the Millennium II simulation (iBovlan-Kolchin et al.1 
120091) . To determine the mass distribution, a power-law ansatz 
was used: 


^=A/o m~ a . 
dm 


(B.4) 


The fit yielded a = 2.0 
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